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ABSTRACT 


An infrared thermographic technique is developed to obtain the transient 
solid surface temperatures surrounding the droplet during vaporization. This 
technique is appealing because it is non-intrusive, detailing the surface 
response to the droplet without affecting the evaporation process. Surface 
recovery can also be monitored using this thermographic method. The 
transient temperature distribution of a Macor solid is detailed. It is found that 
contact temperature is held in the vicinity of the droplet during the majority 
of the droplet’s evaporation until the droplet thickness diminishes greatly, 
where upon the temperature of the solid surrounding the droplet begins to rise. 
The non-dimensional radius of influence of droplet cooling is also detailed. 
The data obtained on the cooling effect induced on aluminum and on Macor 
from previous studies is used in concert with new data obtained on a quartz 
surface to characterize the induced cooling of a hot surface by an evaporating 
droplet. The role of the droplet size and shape is investigated for various 
high and low thermal conductivity surfaces. Droplet evaporation time, surface 


heat transfer coefficient and droplet shape parameter are also examined. 


FOREWORD 


This report describes the research performed during the period July 1988 - 
July 1989 under a joint research program between the Mechanical 
Engineering Department of the University of Maryland and the Center for 
Fire Research of the National Institute of Standards and Technology. The 
research is conducted in the laboratories of the CFR by Mr. Michael Klassen, 
a Graduate Research Assistant of the ME Department under the joint 
supervision of Dr. Marino di Marzo (ME Dept. - UMCP) and of Dr. David 
D. Evans (CFR - NIST). This reports also constitutes the Master Thesis of 
Mr. Klassen which has been completed and defended in the month of 
September 1989. 
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1.0 INTRODUCTION 

This study constitutes a portion of an extensive research program aimed 
at developing an accurate extinguishment model of burning solid fuel surfaces. 
The research deals with the transient cooling induced by the vaporization of 
a single droplet on a hot surface. The majority of the investigation is 
concerned with pure evaporative phenomena (full suppression of nucleate 
boiling), although a limited investigation of the nucleate boiling region is 
presented. 

An infrared thermographic technique is developed to obtain the transient 
solid surface temperatures surrounding the droplet during vaporization. This 
technique is particularly appealing because it is non-intrusive, detailing the 
surface response to the droplet without affecting the evaporation process. 
Surface temperature recovery is also monitored using this thermographic 
method. Furthermore, the extent of droplet influence on the solid can 
accurately measured, as this technique monitors the entire surface of the solid. 

Previous research completed within this program will be utilized for 
comparison. Past data of the cooling effect induced by dropwise evaporation 
on aluminum (di Marzo et al.[1986]) and on Macor® (di Marzo et al. [1987]) 
are used in concert with new data for a quartz surface to characterize the 
effects of material properties. The role of droplet size and shape on droplet 


evaporation time is investigated for the various high and low thermal 


conductivity surfaces. The heat transfer coefficient, evaporation time, and 


shape parameter are evaluated for various surfaces. 


2.0 BACKGROUND 

Liquid sprays are used to enhance the heat transfer from hot surfaces in 
many applications, such as spray evaporators, combustion, fire extinguishment, 
quenching, etc. Many studies have been performed to quantify the 
vaporization process on both single and multi-droplet arrays impacting on hot, 
solid surfaces. These studies generally have dealt with large excess 
temperatures (difference between the surface temperature and liquid saturation 
temperature), though some have spanned the full range of the droplet 
vaporization process, including evaporation, nucleate boiling, film boiling and 
Leidenfrost transition. 

Abu-Zaid and Atreya [1989] studied the transient cooling of a low- 
conductivity ceramic surface by droplet evaporation over a large temperature 
range (75°C to 200°C). The surface was instrumented with embedded 
thermocouples which gave the transient temperature profile along the top 
surface as well as in the depth of the solid. Furthermore, they also reported 
on the recovery time of the surface following the complete vaporization of 
the droplet. 

Seki [1978] measured the transient surface temperature profile of a hot 
wall cooled by a single droplet using a thin-film thermometer. A simple 
solution based on the contact of semi-infinite bodies is presented for the 


change in surface temperature when a cool droplet impacts a hot surface. 


Bonacina [1979] investigated the heat transfer mechanisms during mist 
cooling at relatively low excess temperatures in a high pressure environment. 
This study developed a mathematical model based on dropwise evaporation. 
Rizza [1981] generated a numerical solution for spray evaporation with excess 
temperatures below the Leidenfrost transition temperature and for non- 
isothermal wall conditions. These studies have limited use in the fire 
extinguishment field, though they do establish some common parameters, such 
as the non-dimensional shape factor (B,) used in this report. 

Makino and Michiyoshi have published a number of papers dealing with 
the evaporation of a water droplet on heated surfaces. These papers have 
dealt with the full range of the boiling curve. Furthermore, the materials 
used were high conductivity metals. Among the details discussed in the 
papers are the time-averaged heat flux on the surface and the effect of the 
initial size of the droplet on the evaporation time [1986]. They also described 
a model for the transient heat transfer from the surface to the water droplet 
[1987]. 

This present study is a continuation of an ongoing project to develop an 
accurate droplet cooling model of solid surfaces. Thus, the work completed 
previously is used extensively and built upon. di Marzo and Trehan [1986] 
studied the transient cooling of a hot aluminum plate by water droplets 
evaporation. A computer model to predict the cooling effect induced by a 


single droplet was developed. This work was then used by Evans and di 


Marzo [1986] to derive a model to predict the evaporation of water on a high 
thermal conductivity surface. 

di Marzo et al. [1987] continued the study on a low thermal conductivity 
surface, Macor. Macor is a glass-like solid which can withstand high thermal 
stress. The experimental studies conducted on Macor tabulated evaporation 
times and droplet sizes. A droplet deposition system was developed that 
increased the repeatability of the size and shape of the droplet. This apparatus 


is used in the present study. 


3.0 SCOPE OF WORK 

The solid materials of interest to fire applications exhibit low thermal 
conductivity. As mentioned previously, much of the research in dropwise 
cooling has been conducted on high thermal conductivity surfaces. Thus, an 
objective of the overall project is to extend the previous studies to include the | 
effect of dropwise cooling on low thermal conductivity surfaces. Previous 
work within this research project studied the evaporation of a droplet on 
aluminum and Macor, a low thermal conductivity surface. A large data base 
was compiled for both high and low thermal conductivity surfaces. This data 
base was further supplemented in this study with additional work on a quartz 
surface, another low thermal conductivity material. Attempts to normalize this 
data are introduced in the present study. 

A technique to monitor the behavior of a solid surface during dropwise 
evaporation is developed. Infrared thermography is used as a non-intrusive 
method of detailing the behavior of the solid during the droplet evaporation 
and subsequent recovery to its initial temperature. A greater understanding of 
the physical processes of dropwise cooling is gained by this technique. The 
radius of influence of the droplet is determined for droplets in both the pure 
evaporative regime and in the boiling regime. The results of this study will 
be used as inputs to a multi-droplet model and finally to a solid fuel 


extinguishment model. 


4.0 EXPERIMENTAL SETUP 

4.1 Droplet Deposition 

The apparatus used to deposit the droplet on the surface is shown in 
Figure 4.1. This apparatus was designed to strictly control the droplet release 
height and position of the micro-dispenser with relation to the surface (di 
Marzo et al. [1987]). This degree of accuracy is necessary because these 
factors affect the final configuration (i.e. shape) of the droplet, which is a 
governing parameter affecting the droplet evaporative behavior. 

The solid is heated from below by an electrical heater (Figure 4.2). The 
flux from the heater is controlled by a voltage regulator. The solid is 
surrounded by ceramic fiber insulation to ensure uniform radial temperature 
distribution. Thermocouples are mounted on the solid at both the top and 
bottom surfaces (Figure 4.3) to obtain heat flux data and to assure that the 
solid has achieved a radially uniform temperature distribution. The droplet 
Ae and evaporation are monitored by a video camera with a high magnifying 
lens. The droplet evaporation time is measured with a stopwatch. 

The initial shape of the droplet on the surface and thus the droplet 
behavior is affected by the surface properties and conditions. As every effort 
is made to produce a droplet with a consistent diameter and shape, the surface 
of the solid must be uniform in every respect. To achieve this repeatability, 


the surface is cleaned with 98% ethyl alcohol to remove grease and then with 
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Figure 4-2. Solid Block Modular Enclosure 


* T1, Tz, T3 are ternperature on the bottom 
surtace 


¢ Ts, Ts, Tg and Tg are temperature on the 
top surface 


e T7 is the ambient temperature 


e V; is proportional to the temperature difference 
between the top surface and ambient; and V2 
t8 proportional to the top and bottom surfaces 


Figure 4-3. Thermocouple Location on the Solid 
S 


deionized water. This provides a clean surface and establishes a water 
microlayer prior to droplet deposition. 

A droplet is "gently" placed on the surface of the solid by the dispenser. 
Deionized water is completely degassed through repeated freezing/boiling 
cycles under vacuum (di Marzo et al. [1985]). The water treatment is 
necessary to remove all mineral impurities and dissolved gases from the 
droplet. 

The presence of minerals and gases causes variations in the heat transfer 
characteristics of the droplet. Mineral impurities within the water tend to 
deposit on the solid’s surface causing a degradation of the surface finish. 
Dissolved gases in the water form bubbles within the deposited droplet, thus 
altering the droplet formation (size and shape) and changing the heat transfer 


characteristics of the droplet. 
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4.2 Infrared Apparatus 
The infrared equipment used is the Model 525 Infrared Thermal Imaging 


System made by Inframetrics. Figure 4-4 shows the equipment for taking the 
infrared data. A close-up lens is used to capture the thermographic data of 
the solid within approximately five diameters of the wetted region. A one- 
dimensional temperature profile of the solid surrounding the droplet is found 
using the Line-Scan function of the system (Figure 4-5). This information is 
recorded onto video tape for later retrieval. 

Macor is chosen over quartz as the solid material because of the 
transparency of the quartz solid. This transparency allowed transmission of 
the infrared radiation not only from the surface, but also from within the 
solid. Thus, determining the depth at which the camera is receiving its data 
would be quite difficult. Macor, being an opaque solid, does not present this 
problem. 

The positioning of the infrared camera is important. To achieve the best 
image without picture distortion, the camera should be placed no more than 
45 degrees from the normal of the surface of the solid. Furthermore, the 
distance of the camera from the solid should be maximized to reduce the heat 
input to the infrared equipment from the solid, thus maintaining good quality 
and accuracy of the infrared image. 

A video digitizer is used to acquire images of the temperature distribution 


of the solid from video tape. The images are then used in an image- 
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Figure 4-5 
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processing routine to obtain the necessary information. The video digitizer 
used is Computer-Eyes®, which is a PC-based image acquisition program 
manufactured by Digital Vision. The video tape containing the solid’s 
temperature distribution is paused at certain time intervals, and an image is 
acquired. This information, containing the brightness levels or digital intensity 
values of each picture element (pixel), is then downloaded into a data file. 
The video digitizer recorded the digitized picture into an array of (641 x 201) 
pixel values. Each pixel is assigned a value from 0 to 63 (filling the six bits 
of information) pertaining to its digital intensity, with 0 being black, 63 white 


and the numbers in between being shades of gray. 
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5.0 EXPERIMENTAL RESULTS 


5.1. Infrared Thermography on Macor 
5.1.1. Data Acquisition 


The data received from the infrared equipment is not temperature, but a 
difference in the radiance level between a reference surface of known 
emissivity and temperature and the unknown target surface about which 
information is desired. The reference surface is taken as the solid surface 
at a distance far removed from the droplet perimeter. The general radiometric 
equation is given by: 

Al=(e,*f(T,)+(1-€)*£(T,)} - {e,*£(T,) + (1-€)*£(T,)} G-1) 
where f(T) represents the overall spectral response of the infrared system at 
the temperatures of the various surfaces. This system transfer function is 
unique to each type of infrared system. In the present application, Equation 
5-1 is reduced further, because the emissivity of the target and reference are 
the same, to become: 
Al= ¢€,*{f(T,)-f(T;)} (5-2) 
The background radiance is thus eliminated from the parameters affecting the 
apparent radiance temperature difference. The accuracy of the radiance 
reading by the infrared equipment is further improved when both the target 
and reference surfaces have the same shape or contour and the same 
orientation to the background radiance level. This is true for the present 


application since the reference and target surfaces are taken from the same 


1a) 


solid . Thus, fairly precise and accurate results (within a 20 percent range) 
can be obtained with this equipment and test set-up (Inframetrics Applications 


Manual [1981]). 


5.1.2 Post-test data processing 


a. Theory 


Data received form the infrared thermographer is in the form of a transient 
visual image and is stored on video tape. The goal in the post-test data 
processing is to be able to accurately determine the temperature distribution 
on the solid surface. In order to gain accurate and meaningful results, an 
image processing technique was developed which extracts the solid surface 
temperature distribution data from the image received from the infrared 
thermographic equipment. This temperature data is fit with an exponential 
curve. This is done a number of times during the droplet evaporation, thus 
establishing the transient temperature behavior of the solid. A robust 
technique has been created which can accurately determine the temperature 
profile of the solid for many surface temperatures during the entire droplet 
evaporation process. 

Image processing involves the editing and highlighting of features from 
a digital image, in order to extract useful information. The image received 
from the infrared thermographer (Figure 5-1) contains both the temperature of 
the solid, the droplet, and an isotherm scale. Figure 5-1 is a 641 x 201 


digital representation of the captured image. The only portion of this image 
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that is needed is the temperature of the solid. The temperature of the droplet 
is disregarded, since the water is transparent and it can not be determined at 
what depth the temperature is being scanned. The isotherm scale must be 
eliminated from the picture because these lines tended to interfere with the 
curve fitting of the temperature data. The curve is fitted from the scattered 
image received from the infrared equipment, therefore all other information 
within the image must be removed. 

Image processing begins with the image captured by the video digitizer 
from a paused video tape. In order to extract the necessary information, 
various C language image processing procedures are applied on the data file 
received from the digitizer. Both point and area processes are employed in 
the image processing. A point process works by looking at individual pixels, 
and performing an operation based on that pixel value only. An area process 
works by inspecting the neighborhood surrounding a pixel and performing an 
operation based on the pixel population of that neighborhood [Dawson, 1987]. 

The majority of the unwanted information is discarded by "thresholding". 
Thresholding involves setting a pixel brightness criteria to which all pixels are 
compared (Dawson [1987]). This criterion is determined so that most 
unwanted information is discarded. The use of this technique binarizes the 
data, moving from a data file in which each pixel has a intensity value, to a 
data file which only identifies if the pixels are "on" or "off'. This technique 


generally removes all background portions of the image, leaving the 
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temperature distribution data, any "noise" left in the image, and the isotherm 
grid. 

Erosion is then employed to filter out "noise" in the picture. Since the 
image is acquired from a paused video tape, there tends to be excessive noise 
in the image. Erosion, which is a thinning technique, involves scanning the 
neighborhood surrounding each pixel (those left after thresholding) and 
counting how many pixels are "on" in that neighborhood. The total number 
of "on" pixels in the neighborhood is compared to a previously set value 
which has been shown to remove most extraneous pixels. This procedure 
does not remove the temperature distribution data. This is an interactive 
process as there is no fixed value to filter out all the noise of the image. 
Typically, a comparative value of six "on" elements in a neighborhood 
consisting of 121 elements tends to remove most erroneous points while not 
removing the temperature distribution data. 

The isotherm scale is eliminated from the digitized picture by determining 
the position of each line of the scale and removing those pixels. Each 
isotherm line was an equal distance apart, thus once the top line position was 
obtained, the remaining lines were easily eliminated. However, the position 
of the top line does vary from picture to picture, possibly due to variations 
in the operating conditions when the results were originally recorded onto 


video tape. 
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With all extraneous information removed from the image, the temperature 
distribution of the solid surface is now fitted with a curve (Figure 5-2) using 
a general least-squares fitting routine (Press et.al [1988]). Several different 
types of curves were fit, but it was determined that an exponential fitting was 
generally the most accurate and appropriate. The temperature distribution is 
found have the best fit with an exponential of the form: 

TAS Ave =e (5-3) 
where A and B are constants determined by the fitting routine. The 
orientation of the coordinate system used for these curves is shown in Figure 
5-3. The goodness of fit is determined by the Chi-squared test (Press et 


al.(1988]). 


Figure 5-3 - Coordinate System 
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b. Implementation 


As shown earlier, the data received from the infrared thermographer is in 
terms of the apparent radiance temperature difference (AI) and must be 
converted from this state to actual temperature. Inframetrics included a hand- 
held computer with programs which contained the system transfer function 
mentioned in section 5.1.a. The computer had several programs which could 
convert isothermal difference to absolute temperature. But, the actual transfer 
function could not be extracted from this computer, and efforts to obtain this 
information from the manufacturer were unsuccessful. Thus, AI versus 
temperature curves were generated and curve-fitted for each surface 
temperature using data points obtained from the hand-held computer. These 
are included in Appendix A. The temperature distribution is curve-fitted in 
terms of its screen coordinates, then converted to isothermal units and finally 
converted to temperature using the curves generated from the hand-held 


computer. 


c. Example 


To illustrate the power of this technique, the surface temperature of Macor 
is monitored during the evaporation of a droplet and the subsequent recovery 
of the surface to its original temperature. A 30 tl droplet deposited on the 
Macor solid having a surface temperature of 124°C is selected as an example. 


Figure 5-4 shows the temperature distribution of the solid during the 
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Figure 5-4 - Transient Temperature Distribution of Macor 
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evaporation of the droplet at ten second intervals. The influence of the 
droplet is shown up to 5 radii away from drop center. This information is 
tabulated in Appendix A. Figure 5-5 better shows the extent of the cooling 
process on the Macor solid as evaporation time increases. 

As the transient temperature distribution shows (Figure 5-6), there is an 
initial "shock" to the solid as the droplet first impacts where the surface 
temperature directly surrounding the droplet plunged to approximately 44°C. 
This is expected as a droplet at room temperature (24°C) hits a much hotter 
surface. The droplet quickly then (in approximately 10 seconds) acquires an 
intermediate temperature which is very close to the theoretical contact 
temperature of Seki. Seki [1978] developed a simple solution for determining 
the change in surface temperature when a droplet contacts a hot surface. 
Based on the contacting of two semi-infinite bodies, the contact temperature 
can be given by: 

TeV e aley, 
T. =——-———____ (5-4) 
Yw + Ys 
where y= V pck and T, and T,, are the initial solid and water temperatures. 
For this example, the theoretical contact temperature is 75°C. 

The surface directly surrounding the droplet remains at the intermediate 

temperature of 79°C for approximately 50 percent of the evaporation time and 


then begins to slowly rise as the droplet’s thickness begins to diminish 
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drastically. In contrast, the surface directly surrounding the droplet on a high 
thermal conductivity solid essentially remains at contact temperature for the 
entire evaporation process (di Marzo [1986]). Complete evaporation occurs 
at a temperature close to 10°C above this intermediate temperature. The 
surface recovers in approximately 30 seconds to its initial temperature. As the 
droplet evaporates completely in 101 seconds, the surface recovery is then 
roughly one-third of this time. 

This pattern of solid behavior is consistently seen for all surface 
temperatures in the pure evaporative region. In the boiling region, the solid 
temperature directly surrounding the droplet rises from it’s initial descent to 
approximately saturation temperature where it stays throughout the 
evaporation process. The pattern of surface recovery is maintained. 

The intermediate temperatures for the initial surface temperatures used in 
this study are shown in Figure 5-7. The theoretical contact temperatures are 
found to be in good agreement in the evaporative region, but deviate greatly 
once the boiling region of the curve (approximately 160°C) is reached. This 
is an expected result, since the contact temperatures for the boiling droplets 
should be roughly the saturation temperature of water. This is shown in the 
experimental results. 

Abu-Zaid and Atreya [1989] also monitored the surface temperature of a 
low conductivity ceramic solid by the use of embedded thermocouples. This 


study did not find the drastic initial shock that the infrared thermographic 
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Figure 5-7 - Contact Temperatures of Droplets on Macor 
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technique shows (Figure 5-8). They found the surface under the droplet 
initially falls to approximately contact temperature, then slowly rises to a 
interfacial equilibrium temperature and remains there during the majority of 
the evaporation process. 

The transient change in the non-dimensional radius of influence (r/r,) of 
the droplet on the surface the example is shown in Figure 5-9. The influence 
distance is taken as the distance where the solid’s radiance change is at least 
10% of the maximum change in radiance. The droplet’s influence on the 
solid gradually increases with time until just before complete evaporation, 
where the influence of the droplet jumps dramatically from a value of about 
4.5 to 7. This is due to the exponential increase in the heat flux as the 
droplet thickness diminishes. As the surface recovers towards it’s initial 
temperature, the radius of influence reduces slowly from this maximum value. 
An average radius of influence of 4.1 is found for this particular case. 

The maximum radius of influence of droplets of various volumes at a 
range of differing surface temperatures is presented in Figure 5-10. There is 
quite a bit of scatter within this graph. This can be explained by trend shown 
in Figure 5-9. As shown in this figure, there is a drastic increase in the 
radius of influence as the droplet nears complete evaporation (101 econ 
Thus, if the radius of influence is determined for the droplet at a time 
between approximately 90 and 100 percent of the evaporation time, the 


determined value can be significantly less than the maximum radius of 
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Figure 5-9 - Transient Radius of Influence 
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influence. Since the important parameter for later droplet cooling models will 
be the average radius of influence, this example shows that the radius of 
influence should be determined experimentally at roughly 70 percent of the 
total droplet evaporation time. 

If these values are plotted against the evaporation time of the droplet, a 
pattern seems to emerge for those droplets in the pure evaporative region and 
those in the boiling region. Separate linear regressions are performed for 
droplets in both regions. di Marzo et al. [1987] computed the radius of 
influence for a high thermal conductivity surface and found it increased 
linearly with the evaporation time. A quite different trend was observed on 
the Macor solid. 

The amount of energy required to evaporate a droplet of given volume 
is constant, q= m * h,. This energy is drawn from a volume of the solid 
saapuoding the droplet. The characteristic length of the droplet influence on 
the solid is given by the relation l=Vor. In the evaporative region, the peak 
radius of influence is characterized by a value of about seven. Since the 
evaporation time is relatively long, the evaporation process can draw energy 
from a greater amount of the surrounding solid, thus the radial influence is 
generally greater, but the axial distance influenced by the evaporative process 
is small. In the case of the boiling regime, the short evaporation times lessen 
the length of the droplet influence in the radial direction, but increase the 


axial distance influenced. Thus for the rapid boiling process, the low thermal 


a5) 


conductivity of Macor constrains the surface area affected to a small region 
around the droplet with very intense localized cooling as opposed to the 


evaporation process where more surface area is affected but the temperature 


gradients are less sharp. 
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5.2 Evaporation on high and low thermal conductivity surfaces 


5.2.1 Convective and Radiative heat transfer coefficients 

The overall heat transfer coefficient is evaluated by measuring the heat 
flux through the semi-infinite solid and the solid-surface-to-ambient-air 
temperature difference. The details of these experimental measurements are 
described in previous work (di Marzo et al. [1986] and [1987]) and later in 
this paper. 

The overall heat transfer coefficient is easily found for the Macor and the 
aluminum solids since the temperature of the solid is known at both the top 
and bottom surfaces. This data are tabulated in Appendix B. 

The calculation of the convective heat transfer coefficient for the quartz 
solid is not a trivial task. Since this is a transparent medium, the internal 
reflections make this calculation much more difficult (Viskanta and Anderson 
[1975]) and was not attempted at this time. 

The convective heat transfer coefficient is deduced from the overall heat 
transfer coefficient data. The convective heat transfer coefficient is found by 
realizing that the heat is removed from the solid surface by a combination 
of convection and radiation according to the following relationship: 

He(@lee be). Clee Ol eee) (5-5) 
where the emissivity of the surface are e= 0.17 and e= 0.84 for aluminum 
and Macor respectively. This equation can be rearranged if a radiative heat 


transfer coefficient is defined: 
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h= eo(T,” +T,’)(T, + T,) (5-6) 


The convective heat transfer coefficient can now be obtained from equations 
(5-5) and (5-6). The results are given in Appendix B. Theoretical values of 
the heat transfer coefficient are found by considering the exposed surface as 
a hot disk facing upward. This can be done since the sides of the solid are 
well insulated and hence the heat loss from the sides is negligible. A 
correlation for the convective heat transfer coefficient in natural convection 


from a hot disk facing upward is proposed by McAdams [1957]: 


T, 
h.= 1.32 (———)* ay) 

S/P 
where S/P is the ratio of the disk surface area to its perimeter. This 
correlation gave consistently low theoretical values compared to the 
experimental results (di Marzo et. al. [1989]). Further investigation showed 
that there is a wide range of experimental correlations for the heat transfer 
coefficient of a hot disk facing upward. Goldstein et al. [1973] compiled 
many of these empirical results. Figure 5-11 plots five different empirical 
relations presented in Goldstein et al. [1973] and the convective heat transfer 


coefficients computed from this study. This figure shows that the coefficients 
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Figure 5-11 - Experimental Convective Heat Transfer Coefficients 
versus Emperical Values 
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calculated in this study are well bounded by these published empirical 


relations. In particular, the relation: 


h.=.59(Ra)’"*(k/L) (5-8) 


where L again is the ratio of the surface area to the perimeter fits this study’s 
experimental data well. 

The fact that the experimental values found in this study are well bounded 
by the various empirical relations establishes more confidence in the validity 
of these results. The heat transfer coefficients were determined from 
experimental results compiled over several years by different people. Since 
this data is well bounded, these findings show that the procedure is reliable 
and robust. 

3.2.2 Evaporation Time 

The results discussed here are a compilation of the aluminum (di Marzo 
et al. [1986]) and Macor (di Marzo et al. [1987]) data found previously, along 
with additional data found for quartz surface. The droplet sizes used ranged 
from 10 to 50 wl. An average of ten droplets are used at each different 
volume to obtain the final value. The quartz data are presented separately, 
and then compared with the aluminum and Macor data in order to make 
qualitative statements about the cooling of a hot surface by dropwise 


evaporation. 
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Figure 5-12 shows the evaporation times for the various surface 
temperatures and droplet volumes on the quartz solid. The shape of the curve 
is similar to those found in the earlier work on Macor and aluminum. The 
actual average evaporative times for each surface temperature and volume are 
given in Appendix B. The evaporation times are compared for various 
surfaces in Figure 5-13. If the evaporation time is graphed on a natural log 
scale, a fairly straight line can be obtained. Figure 5-14 shows this scaling 
for the quartz data, while a compilation for the three surfaces at a volume of 
30 tl is shown in Figure 5-15. 

If logarithmic regression is performed on the data from each surface, the 
evaporation time can be modelled in the following form: 

T=1A.e. 5% (5-9) 
where the constants for each surface for various droplet volumes are found in 


Table 5-1. The constants for the aluminum solid are from di Marzo [1986]. 
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Table 5-1 - Coefficients of Fitted Evaporation Time 
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Figure 5-12 - Droplet Evaporation Times on Quartz Solid 
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Figure 5-13- Droplet Evaporation Times for Different Solids 
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Figure 5-14 - Natural Log of Droplet Evaporation Times 
for Quartz Solid 
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Figure 5-15 - Natural Log of Droplet Evaporation Time 
all Solids 
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Figure 5-16 shows the average shape parameter (8) found for each droplet 
volume at each initial surface temperature for the quartz solid. The values for 
this graph are compiled in Appendix B. Great care was taken to place a 
droplet of consistent shape and size on the solid. If the 8) values for each 
volume are averaged at each surface temperature (Figure 5-17), the shape 
parameter is shown to average 1.53 with a variance of approximately three 
percent. 

The shape parameter is greater for the quartz surface averaging 1.53 as 
compared to 1.42 that di Marzo et al [1987] found for the Macor solid 
(Figure 5-18). This means that the size of the wetted region on the quartz 
was consistently larger, thus giving a larger wetted area of contact to the 
solid. It is interesting to note that at the corresponding surface temperature, 
evaporation time is longer on the quartz surface than on the Macor surface. 
Both surfaces have similar thermal properties (Appendix B), thus this finding 
causes one to question the relative influence of B, on evaporation time. 
Bonacina [1979] found that the shape parameter was very influential in the 
evaporation time of a droplet, on the order of three. The results of this study 
seem to be at variance with this assertion. This point is supported further by 
Figure 5-19. This figure contains evaporation times for all cases, including 
those in which f, was not strictly controlled. B, ranges from 1.43 to 1.7 for 


the droplets included in this curve. Though there is some fluctuation, the 
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Figure 5-16 - Droplet Shape Parameters for Quartz Solid 
at each Droplet Volume 
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Figure 5-17 - Average Droplet Shape Parameter for Quartz Solid 
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Figure 5-18 - Comparison of Droplet Shape Paramenter 
on Different Solids 
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Figure 5-19 - Droplet Evaporation Times for Large Range of 8 
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general trend of the curve is similar to that found when control of B, is 
maintained. 

5.2.3 Evaporation Time in Boiling Regime 

A limited study of the evaporation of droplets in the boiling regime is 
undertaken on both the quartz and the Macor solids. The trend of the 
evaporation time curves found for the pure evaporative region are extended 
by this data. 

Figure 5-20 shows the evaporation times on both quartz and Macor. 
Appendix B contains these values. Control of the shape of the droplets is 
quite difficult because the droplet tends to spread after impact. Presumably 


this is due to lessening of the effect of surface tension at these temperatures. 


5.2.4 Normalization of Evaporation Time 


Droplet evaporation is governed by numerous parameters which affect the 
menoeiiccid interactions. Among these parameters are the thermal conductivity | 
of the solid, the thermal diffusivity of the solid, and the shape parameter. 

As shown by the infrared thermography, the temperature of the solid 
surrounding the droplet tends to remain constant for the majority of the 
droplet evaporation time. This so-called contact temperature is an important 
parameter when considering dropwise evaporation. As introduced before, the 
contact-temperature can be calculated fairly accurately in the pure evaporative 


region using equation 5-4. If the evaporation times for all the surfaces are 
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Figure 5-20 - Droplet Evaporation Times in Boiling Regime 
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plotted against their contact temperatures (Figure 5-21), the curves tend to 
collapse along the horizontal axis. The collapse seen is expected, as the 
contact temperature should remain around 100 °C as the droplet enters the 
nucleate boiling regime. This trend shows that the contact temperature is an 


important parameter in reducing the evaporation process for all surfaces. 
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Figure 5-21 - Normalization of Droplet Evaporation Times 
by Contact Temperature 
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6.0 SUMMARY AND CONCLUSIONS 

6.1 Summary 

An infrared thermographic technique has been developed to obtain the 
transient solid surface temperatures surrounding the droplet during 
vaporization. The surface of the solid surrounding the droplet can be 
monitored, without affecting the evaporation process. Surface recovery to its 
initial state can also be monitored. The transient temperature distribution of 
a Macor solid is detailed. The non-dimensional radius of influence of droplet 
cooling is also found. 

The data obtained on the cooling effect induced on aluminum and on 
Macor from previous studies is used in concert with new data obtained on a 
quartz surface to characterize the induced cooling of a hot surface by an 
evaporating droplet. In particular, the role of droplet size and shape is 
investigated for the various high and low thermal conductivity surfaces. A 


limited study of droplet vaporization in the boiling regime is also presented. 


6.2 Conclusions 

The following conclusions can be drawn from this study: 

The infrared thermographic technique improves the physical understanding 
of the processes which take place during dropwise evaporation. In particular, 
it is found that solid surface drops to an intermediate temperature, which is 


well approximated by Seki’s theoretical contact temperature. This contact 


fo 


temperature is held in the vicinity of the droplet edge during the majority of 
the droplet evaporation until the droplet thickness diminishes greatly, where 
upon the temperature of the solid surrounding the droplet begins to nse. The 
region affected by the droplet evaporative cooling is characterized with the 
non-dimensional radius of influence which is determined using the infrared 
thermographic technique. It is found that the radius of influence increases 
dramatically near the end of the evaporation process. When determined over 
a range of surface temperatures, the peak radius of influence is found to be 
approximately seven wetted region radii for the pure evaporative region. In 
the boiling region, the radius of influence varies within the range of 2 to 5 
wetted region radii. 

The additional study of droplet evaporation on a quartz solid has improved 
the data base started with the previous studies on aluminum and Macor. The 
average droplet shape parameter is greater for the quartz (1.53) as compared 
to 1.42 found previously on the Macor. However, it is found that for the 
corresponding surface temperature, the evaporation time on the quartz surface 
is longer than that found on the Macor surface, though they have similar 
material properties. This finding raises questions on the relative importance 
of the shape parameter on droplet evaporation time. The study of droplet 
evaporation in the boiling regime showed the general trend found in the pure 


evaporative region is extended. 
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The results of this study will be used as inputs for a multi-drop model of 
the transient cooling of a hot surface by dropwise evaporation. Subsequent 
studies will look into dropwise cooling of a surface heated by radiation, thus 


representing a more realistic fire environment. 
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Image Processing Programs 
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Program Lineout a 


/* INCLUDES 
include <stdio.h> 
include <conio.h> 
include <math.h> 
include <stdlib.h> 
include <dos.h> 
include <graphics.h> 
include <string.h> 
include "nrutil.h" 
include "nr.h" 
include <mouse.inc> 
include <gmouscur.inc> 


HHHtHHHAAHHAH +H 


# include "“arry.h" 


define ROW 201 
define COL 320 
define LEFT 0 
define RIGHT 1 
define NTERM 2 


te te te Ht H+ 


/* LOCAL PROTOTYPES 
unsigned char ** read_image(char *); 
void graphicScreen(void); 
void Re_Order(unsigned char **); 
void get_range(void); 
void check_hood(int,int,int); 
void Isf(float*,float*,int, float *, int, FILE *); 
void free_vector(float*,int,int); 
float *vector(int,int); 
void Boarder(void); 
void Draw(float *, float *, int, int,int); 
/* GLOBAL VARIABLES 


unsigned char **a; 

float *x,*y,*aa; 

int co,ro,cl,rl; 

int d1,d2,e1,e2,dropcntr,offset; 
FILE *Stats; 


void main() 
unsigned char **v; 
char in_file[25]; 


char pre [30],pre2 [30]; 
register int i, j; 
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va 


*) 


rf 


fs 


int Tr, Max,countr,flag,Nr,droprad; 
FILE *Outfile; 


/* VARIABLES ra} 
/*v = vector read in from image file. lf 
Hb Must be reordered since PC-Eyes stores image */ 
/* in vertical columns, but C reads across rows.*/ 

/*a= reordered matrix */ 


a= (unsigned char **)matrx(ROW, 2*COL,sizeof(unsigned char *),sizeof(unsigned char)); 


x=Vvector(1,1000); 
y=vector(1,1000); 
printf("Enter the image filename -->"); 


gets(in_file); 

strcpy(pre,in_file); 

strcat(pre,".dat"); 

Outfile= fopen(pre,"w+"); 

strcpy(pre2,in_file); 

strcat(pre2,"st.dat"); 
Stats=fopen(pre2,"w+"); 


fprintf(Stats,'\n",in_file); 
Vv = read_image(in_file); 


Re_Order(v); 


gotoxy(400,10); 

printf("Threshold? -->"); 

scanf("%d", &Tr); 
fprintf(Stats,"\nThreshold= %d\n",Tr); 
printf("\nMax? -->"); 

scanf("%d" ,& Max); 

fprintf(Stats,"Max= %d\n ",Max); 
printf("Number of rows to first line -->"); 
scanf("%d" ,& Nr); 

fprintf(Stats,"Number of rows to first line= %d\n",Nr); 
printf("Offset--> ",offset); 
scanf("%d" ,& offset); 

fprintf(Stats,"Offset= %d\n" offset); 


Thresholding and neighborhood thinning */ 


for(i = 0; i < ROW; i++) 
for(j = 0; j < 2*COL; j++) 
if(i<19)a[i) [jJ=0; 
else if (2>195)a[i)[j]=0; 
else if(a{i)[j] >= Tr) 
a{iJ{j) = 1; 
else 


aliJ) = 0; 
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for(i = Nr; i <ROW; i=i+16) 
for(j = 0; } < 2*COL; j++) 


{ 
a[i}[j)=0; 
ali+1]j]=0; 
} 


for(i = 0; i < ROW; i++) 
for(j = 0; j < 2*COL; j++) 
if (afiJJ==1) 
check_hood(i,j,Max); 


graphicScreen(); 
Boarder(); 


for(i = 0; 1 < ROW; i++) 
for(j = 0; } < 2*COL; j++) 
if(a{iJQ] == 1) 
putpixel(j, 1, 6); 


jf Select range of Curve =] 


get_range(); 
fprintf(Stats,"The data was taken from the range %d %d\n",co,c1); 
fprintf(Stats,"The droplet is between %d %d",d1,d2); 


/* Selecting values from which Polynomial is fitted “fj 
countr=0; 
droprad=abs(dropcntr-d1); 


for(i=ro; i<=rl; i++) 
for(j=co; j<= C1; j++) 

if(afiJUJ==1) { 

countr+=1; 

x[countr]=(float) j; 

y(countr]=(float) 1; 

y(countr]=log(y[countr]); 

if(dropcntr>=x[countr]) 


x{countr]=dropcntr+offset-x[countr]; 
flag=0; 
) 


else if(dropcntr<x[countr]) 


x{countr]=x[countr]-dropentr-offset; 
flag=1; 
) 


} 
forG=1; j<= count; j++)fprintf(Outfile,"x= %4.5f  y=%4.5f\n",x[j],-y[j]); 
fprintf(Outfile,"droprad= %d\n ",droprad); 
/*  fprintf(Outfile,"countr= %d ",countr);*/ 


free_matrx((void**)a,ROW); 


[EE EEEEEEEREL ER ELE. EEE EERE S EERE EEE REE RS EEE RW Se ek ee 


{= Calling Least Squares Fit Rouune */ 
=vector(1,NTERM); 


Isf(x,y,countr, aa, NTERM, Stats); 


[* Calling Plotting Routine ah 
y=vector(1,countr+ 10); 
for(i = 1; i <= countr; i++) 
{ 
yli] = 0.; 


y{i] += exp(aa[1])*exp(-aa[2]/x[1]); 


} 
/* forgG=1; j< counu,; j++) fprintf(Outfile,"x= %4.5f y=%4.5f\n".x{jJ,yUj));*/ 


Draw(x, y, count, 1,flag); 


fclose(Outfile); 
fclose(Stats); 


free_vector(aa, 1, NTERM); 


free_vector(x,1,1000); 
free_vector(y,1,countr+1); 


outtextxy(400, 300, "HIT RETURN"); 
getch(); 
closegraph(); 


void get_range(void) 
{ 
resetRec *theMouse; 
locRec *L; 
int count = 0; 
char cbuf[50]; 


/* Set up for a run. = 
initGCurs(); /* Initialize the cursor images. = 
theMouse = mReset(); /* Reset the mouse. og 


if(theMouse->exists) 


{ 
outtextxy(480,10, "Cursor Position”); 
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/* Show the default cursor. 


setcolor(10); 
mShow(); 

/* Set the range for the cursor. 
mColRange(0,2*COL); mRowRange(0,ROW); 

/* Select the hand cursor. 
mGraphCursor(hand.hotX, hand.hotY, (unsigned)(_DS), 

(unsigned)(hand.image)); 

gotoxy(60,1); 

/* Wait for a click of the right button. 


do { 
L = mReleased(LEFT); 
if(L->opCount > Q) 
{ 


if(f{mod(count,4.) == 0.) 
{ 

co = L->column; 
ro = L->row; 
mHide(); 
circle(co,ro,5); 
sprintf(cbuf,” 48 
Outtextxy(472, 16, cbuf); 
sprintf(cbuf, "x0 = %3d yO = %3d\n", co, ro); 
Outtextxy(472, 16, cbuf); 
mShow(); 


} 
if(fmod(count,4.) == 1.) 
{ 


cl = L->column; 

tl = L->row; 

mHide(); 

circle(cl,r1,5); 

sprintf(cbuf,” "i 
outtextxy(472, 24, cbuf); 

sprintf(cbuf, "x1 = %3d yl = %3d", cl, rl); 
outtextxy(472, 24, cbuf); 

mShow(); 


} 
if(fmod(count,4.)== 2.) 
{ 


d1l=L->column; 
el=L->row; 
mHide(); 
circle(d1,e1,3); 
mShow(); 


} 
if(fmod(count,4.) == 3.) 
{ 
d2=L->column; 
e2=L->row; 


mHide(); 
circle(d2,e2,3); 
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mShow(); 
} 


count += 1; 


} 


L = mReleased(RIGHT); 
}while(L->opCount == 0); 
dropcnir=(d1+d2)/2; 
fprintf(Stats,"droplet center= %d",dropcntr); 
theMouse = mResetQ); 


} 


/* Reordering vector coming in from PC-Eyes 


void Re_Order(unsigned char **v) 


{ 


iInGinys 


for(k = 0; k < 2; k++) 
for(j = 0; j} < COL; j++) 
for(i = 0; i < ROW; i++) 
a{iJ[j+ COL*k) = v(k][j*ROW+i]; 
free_matrx((void **)v, 2); 


/* Thinning routine 


void check_hood(int m, int n, int B) 
{ 
int iJ; 
int count,d=5,e=5,e1=5,di=5; 
count=0; 
if (m<e) e=m; 
else if(n<d) d=n; 
else if ((m+el)>ROW-1) el=(ROW-1)-m; 
else if ((n+d1)>2*COL-1) dl=(2*COL-1)-n; 
/* — pmintf("e,d,el,dl= %d.%d,%d,%d\n" ,e,d,e1,d1);*/ 
for (i= m-e; i<= m+el; i++) 
for(j=n-d; j<=n+d1; j++) { 
if(a{i]J==1) 


count+=1; 


63 


vol 


if (count<B) a[m][n] = 0; 


/* Set up the graphics screen. dt i 
void graphicScreen(void) 


{ 
int driver = DETECT, mode; 


initgraph(&driver, &mode, "d:\\c"); 
if(graphresult() == grOk) 


setcolor(4); 
} 
} 
i Be oc eeccceswccace ec oe ne eee < 
[= Boarder 4/27/89 
I fee es ei es che wi cbes wc wits ane creo ene eS en ee ee ee en eo a | 


void Boarder(void) 
| /* Draw a boarder around the graph. 
rectangle(0, 0, getmaxx(), getmaxy()); 
/* Draw the x and y axis. 


line(getmaxx()/2, 0, getmaxx()/2, getmaxy()); 
line(O, getmaxy()/2, getmaxx(), getmaxy()/2); 


{Ces oe ie REN Ale Rear 5 oN a i ec per eter nt */ 
tbe Draw 4/28/89 

(eile he EERIE) Sh DT SPRAIN oy nce a CE ri nc */ 
/* 


void Draw(float *x, float *y, int n, int col,int flag) 


int i; 
char buf[25]; 


Hp 


S 
= 


*/ 


Wl 


setcolor(col+1); 
sprintf(buf, "Data File #%2d", col); 
outtextxy(50, 50 + 10*col, buf); 
for(i=1;i<=n;i++) 
if(flag==0) 
x(iJ=dropcentr+offset-x [i]; 
else if(flag==1) 
x(iJ=dropentr+offset+x [i]; 


/* Plot the real part of y. 
for(i = 0; i <= n; 1 ++) 


{ 


iG == 0) 
moveto(x{1], y[1]); 
else 
lineto(x[i], y[i]); 
} 


/* Scale to the screen coordinates. 
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*/ 


/*Subroutine to read in image file */ 


#define ROW 201 
#define COL 320 
#include <stdio.h> 
#include <fcntl.h> 
#include <sys\stat.h> 
#include <io.h> 
#include <string.h> 
#include <process.h> 
#include "nruul.h" 
#include "nr.h" 


#include "arry.h" 


unsigned char **read_image( char *in_file ) 


int fp_in; 

register int 1; 

char pre(30]; 
unsigned char **a; 


strcpy(pre,"a:\\"); 
surcat(pre,in_file); 
strcat(pre,".raw"); 
printf('\n%s\n" pre); 


fp_in=open(pre,O_RDONLYIO_BINARY); 
if(fp_in <0) 
{ 
perror("\nCan’t Open In_File"); 
exit(0); 
} 
a = (unsigned char **)matrx(2, ROW*COL, sizeof(unsigned char*), sizeof(unsigned char)) ; 


for(i = 0; 1 < 2; i++) 
read(fp_in, &a[iJ[0], COL*ROW); 


close(fp_in); 
retum(a); 


JAZ 


ip: Pon ed 


PARP ERE EEE ELAR EERE SS EEL E SEE EE CEE ERE 


tbs Least Squares Fit sal | 

he 06/09/89 dt 

ORS ee ee ee Ke SE eee eS CEH f 
#include <stdio.h> 

#include "nr.h" 

#include “nruul.h" 

#include <math.h> 

#define ROW 201 


i Local Prototypes | 
float *vector(int,int); 

float **matrix(int,int,int,int); 

int *ivector(int,int); 

void free_vector(float*,int,int); 

void free_matrix(float**,int,int,int,int); 

void free_ivector(int*,int,int); 


/*void Ifit(float*,float*,float* ,int,float*,int,int*,int,float**,float**,void**);*/ 


7% Basis Funcuon Subroutine 


void funcs( x, afunc, ma) 


float x,*afunc; 


int ma; 

( 
int 1; 
afunc[1]=1.0; 


for(i=2;i<=ma;i++) afunc[i]=-1/x; 


} 

I ps Dok mesin ns da RED ae DORE a Pen RE REECE RGAE EER OEE eae Sot aeduneaan od) 

fii Driving Program 

is Variables 

/* sig=standard deviations 

ibe chisq=Chi squared test 

/* lista= coefficients to be evaluated el, 
/* covar= covariance matrix “a 
/* Ifit= fitting routine 

| . - pa ada UGS Ss SnmaaneaTESSGEne THESE SEMEEScuauaearsasancosntenseonnenases */ 


void Isf(float *x, float *y, int NPT, float *aa, int NTERM, FILE *Stats) 
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be main */ 
float *sig,**covar,chisq; 
int *lista,ij,mfit; 


lista=ivector(1,NTERM); 
sig=vector(1,NPT); 
covar=matrix(1,NTERM,1,NTERM); 


for (i=1;1<=NPT;i++) 
sig[iJ=1; 


for (i=1;i<=NPT;i++) printf("sig=%2.4f" sig[i]);*/ 
mfit=NTERM; 
for(i=1;i<=mfit;i++) lista[i]=1; 


Calling Fitting Routine fl | 
Ifit(x,y,sig, NPT,aa,NTERM, lista,mfit,covar,&chisq,funcs); 


Printing out Results nef 


fprintf(Stats,"\n%11s %21s\n","parameter","uncertainty"); 
for (i=1;1<=NTERM;i++) 
fprintf(Stats, "aa[%d]=%8.6f %12.6f\n",i,aa[i},sqrt(covar(i][i]})); 
fprintf(Stats,"Chi-squared = %12f\n"chisq); 
fprincf(Stats,"Full covariance matrix\n"); 
for(i=];i<=NTERM;i++) { 
for(j=1;j<=NTERM;j++)fprintf(Stats,"%12f" ,covar[i)[j]); 
} 


free_matrix (covar,1,NTERM,1,NTERM); 
free_vector(sig,1,NPT); 
free_vector(y,1,ROW); 
free_ivector(lista,1,NTERM); 
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the Lfit.c 
Subroutune from Numerical Recipes (Press et. al. [1988]) to perform Least Squares 
|S Biase i 


static float sqrarg; 
#define SQR(a) (sqrarg=(a),sqrarg*sqrarg) 


Local Prototypes =; 


void gaussj(float**,int,float**,int); 
void covsrt(float**,int,int* int); 

void free_vector(float*,int,int); 

void free_matnx(float**,int,int,int,int); 
float *vector(int,int); 

float **matrix(int,int,int,int); 

void *(funcs)(float,float *,int); 

void nrerror(char*); 


void Ifit(x,y,sig,ndata,a,ma,lista,mfit,covar,chisq,funcs) 

int ndata,ma,|lista{],mfit; 

float x{J,y[),siz(J,a{),**covar, *chisq; 

void (*funcs)(); /* ANSI: void (*funcs)(float,float *,int);*/ 


int k,kk,j,ihit,1; 

float ym,wt,sum,sig2i,**beta,*afunc; 

void gaussj(),covsrt(),nrerror(),free_matnx(),free_vector(); 
float **matrix(),*vector(); 


beta=matnx(1,ma,1,1); 
afunc=vector(1,ma); 
kk=mfit+1; 
for §=1;j)<=ma;j++) { 
ihit=0; 
for (k=1;k<=mfit;k++) 
if (lista{k] == j) ihit++; 
if (ihit == 0) 
lista{kk++]=J; 
else if (ihit > 1) nrerror("Bad LISTA permutation in LFIT-1"); 


} 
if (kk != (ma+1)) nrerror("Bad LISTA permutation in LFIT-2"); 
for (j=1;j<=mfit;j++) { 

for (k=1;k<=mfit;k++) covar[j][k]=0.0; 

beta[j]{1}=0.0; 


for (i=1;i<=ndata;i++) { 
(*funcs)(x[i],afunc,ma); 
ym=y(i); 
if (mfit < ma) { 
for (j=(mfit+1);j<=ma;j++) 
ym -= a{lista{j])*afunc[lista[j])}; 
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} 


} 
sig2i=1.0/SQR(sig[i]); 
for (j=1;j<=mfit;j++) { 
wt=afunc[lista[j]]*sig2i; 
for (k=1;k<=j;k++) 
covar[j][k] += wt*afunc[lista[k]]; 
beta[j][1] += ym*wt; 
} 


} 
if (mfit > 1) 

for (j=2;j<=mfit;j++) 

for (k=1;k<=j-1;k++) 
covar[k][j]=covar[j)[k]; 

gaussj(covar,mfit,beta,1); 
for (j=1;j<=mfit;j++) a[lista[{j]}=beta[j][1]; 
*chisq=0.0; 
for (i=1;1<=ndata;i++) { 

(*funcs)(x[i],afunc,ma); 

for (sum=0.0,j=1;j<=ma;j++) sum += a[j]*afunc[j]; 

*chisq += SQR((y[i]-sum)/sig[i}); 
} 
covsrt(covar,ma,lista,mfit); 
free_vector(afunc,1,ma); 
free_matrix(beta,1,ma,1,1); 


#undef SQR 


ae 
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/* Nrutil.c 
Subroutine from Numencal Recipes to dynamically allocate vectors and Matrices */ 


#include <stdio.h> 
#include <stdlib.h> 
#include <alloc.h> 


void nrerror(char*); 


void nrerror(error_text) 
char error_text[]; 


void exit(); 


fprintf(stderr,"Numerical Recipes run-time error...\n"); 
fprintf(stderr," %s\n" ,error_text); 

fprintf(stderr,"...now exiting to system...n"); 

exit(1); 


float *vector(nl,nh) 
int nl,nh; 


{ 


float *v; 


v=(float *)malloc((unsigned) (nh-nl+1)*sizeof(float)); 
if ('v) nrerror("allocation failure in vector()"); 


return v-nl; 

} 

int *ivector(nl,nh) 

int nl,nh; 

{ 
ine ve 
v=(int *)malloc((unsigned) (nh-nl+1)*sizeof(int)); 
if (!v) nrerror("allocation failure in ivectorQ"); 
retum v-nl; 

} 

double *dvector(nl,nh) 

int nl,nh; 

{ 
double *v; 


v=(double *)malloc((unsigned) (nh-nl+1)*sizeof(double)); 
if (!v) nrerror("allocation failure in dvector()"); 
retum v-nl; 


float **matix(nrl,nrh,ncl,nch) 
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int nrl,nrh,ncl,nch; 


int 1; 
float **m; 


m=(float **) malloc((unsigned) (nrh-nrl+1)*sizeof(float*)); 
if (!m) nrerror(“allocation failure 1 in matnx()"); 
m -= nr; 


for(i=nrl;i<=nrh;i++) { 
m{iJ=(float *) malloc((unsigned) (nch-ncl+1)*sizeof(float)); 
if (!m[i]) nrerror("allocation failure 2 in matrix()"); 
m[i] -= ncl; 

} 


retum m; 


} 


double **dmatnx(nrl,nrh,ncl,nch) 
int nrl,nrh,ncl,nch; 
tee 

int 1; 

double **m; 


m=(double **) malloc((unsigned) (nrh-nrl+1)*sizeof(double*)); 
if (!m) nrerror("allocation failure 1 in dmatnx()"); 
m -= nr; 


for(i=nrl;i<=nrh;i++) { 
m{i]=(double *) malloc((unsigned) (nch-ncl+1)*sizeof(double)); 
if (!m[i]) nrerror("allocauon failure 2 in dmatrix()"); 
m{i] -= ncl; 

} 


retum m; 


} 


int **imatrix(nrl,nrh,ncl,nch) 
int nrl,nrh,ncl,nch; 


{ 

Hifai.s ml 

m=(int **)malloc((unsigned) (nrh-nrl+1)*sizeof(int*)); 

if (!m) nrerror(“allocation failure 1 in imatrix()"); 

m -= nrl; 

for(i=nrl;i<=nrh;i++) { 
m{i]=(int *)malloc((unsigned) (nch-ncl+1)*sizeof(int)); 
if (!m[i)) nrerror("allocation failure 2 in imatrix()"); 
m{i] -= ncl; 


retum m; 


float **submatrix(a,oldrl,oldrh,oldcl,oldch,newrl,newcl) 
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float **a; 
int oldrl,oldrh,oldcl,oldch,newrl,newcl; 


int i,j; 
float **m; 


m=(float **) malloc((unsigned) (oldrh-oldrl+1)*sizeof(float*)); 
if (!m) nrerror("allocation failure in submatrix()"); 
m -= newrl; 


for(i=oldrl j=newrl;i<=oldrh;i++,j++) m[j]=a[i]+oldcl-newcl; 


retum m; 


void free_vector(v,nl,nh) 
float *v; 
int nl,nh; 
{ 

free((char*) (v+nl)); 
} 


void free_ivector(v,ni,nh) 
int *v,nl,nh; 
{ 


free((char*) (v+nl)); 


void free_dvector(v,nl,nh) 
double *v; 
int nl,nh; 


free((char*) (v+nl)); 


void free_matrix(m,nrl,nrh,ncl,nch) 
float **m; 

int nrl,nrh,ncl,nch; 

{ 


int 1; 
for(i=nrh;i>=nrl;i--) free((char*) (m[i]+ncl)); 
free((char*) (m+nrl)); 


} 


void free_dmatrix(m,nrl,nrh,ncl,nch) 
double **m; 
int nrl,nrh,ncl,nch; 


{ 


int 1; 
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for(i=nrh;i>=nrl;i--) free((char*) (m[iJ]+ncl)); 
free((char*) (m+nrl)); 


void free_imatnx(m,nrl,nrh,ncl,nch) 
int **m; 
int nrl,nrh,ncl,nch; 


{ 


int 1; 


for(i=nrh;i>=nrl;i--) free((char*) (m[i]+ncl)); 
free((char*) (m+nr]l)); 


void free_submatrix(b,nrl,nrh,ncl,nch) 
float **b; 
int nrl,nrh,ncl,nch; 


free((char*) (b+nrl)); 


float **convert_matrix(a,nrl,nrh,ncl,nch) 
float *a; 
int nri,nrh,ncl,nch; 


int 1,j,nrow,ncol; 
float **m; 


nrow=nrh-nrl+1; 

ncol=nch-ncl+1; 

m = (float **) malloc((unsigned) (nrow)*sizeof(float*)); 
if (!m) nrerror("allocation failure in convert_matrix()"); 
m= nri: 

for(i=0,j=nrl;i<=nrow-1;1++,j++) m{j]=a+ncol*i-ncl; 
return m; 


void free_convert_matrix(b,nrl,nrh,ncl,nch) 
float **b; 
int nrl,nrh,ncl,nch; 

free((char*) (b+nrl)); 


AZ 
AZ 
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Procedure for Obtaining Radius of Influence: 


i Acquire images from video tape using digitizer. It is necessary to 
determine the droplet diameter during the majority of the evaporation process, 
SO acquire an image fairly early (about 30 seconds) into the process. For the 
peak radius of influence, acquire the image as close to the final evaporation 


as possible. 


2. Use image processing program Lineout.exe to fit exponential curve to 
temperature profile. Use this program to also find the droplet diameter in 
screen coordinates from the early image. The maximum radiance drop is also 


found from this image. 


3. Using the coefficients and other data from least squares fitting, such as 
drop radius and drop center, compute the r and z values of curve of the 
images. These will be in screen coordinates. Convert z_ values to 
isotherm units (those of infrared equipment). Divide radial values by radius 


of droplet (this is radius of influence). 


4. If radius of influence is taken as 10% of maximum radiance 
difference of the image acquired at 30 seconds, find where axial value is 


1/10th of maximum radiance difference and note radius of influence. 


Tp 


Procedure for Converting Radiance Difference to Temperature: 

1. Convert axial positions from screen coordinates to isothermal. 

2. From hand-held computer, generate curves of isothermal difference vs. 
temperature. This must be done for each surface temperature. Fit curves to 


these lines. 


3. Input isothermal values found in step 1 into curves. 
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Surface temperature (deg C) 


200 


100 


-100 


y = 100.85 + 9.1115x - 1.5749x%2 - 0.34364x*3 - 3.3205e-2x%*4 


Surface Temperature 101 


Ti 


Surface Temperature (deg C) 


200 


100 


-100 
-12 


y = 123.77 + 8.5577x - 0.89895x4’2 - 0.15323x*3 - 1.2290e-2x*4 


Surface Temperature 124 
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Surface Temperature (deg C) 


160 


140 


120 


100 


80 


60 


40 


20 


y = 142.89 + 9.6490x + 8.3100e-2x%2 + 2.0785e-2x*3 - 6.9930e-4x*4 


Surface Temperature 143 
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Surface Temperature (deg C) 


160 


140 


120 


100 


80 


60 


y = 152.01 + 8.6236x - 0.35015x*%2 - 5.1807e-2x*3 - 3.9336e-3x%4 


Surface Temperature 152 
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Surface Temperature (deg C) 


180 


160 


140 


120 


100 


80 


60 


40 


y = 160.02 + 9.0011x - 2.6515e-3x’2 + 4.5260e-3x*3 - 8.4499e-4x*4 


Surface Temperature 160 
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Surface Temperature (deg C) 


200 


100 


y = 177.92 + 19.438x - 1.4814x*2 - 0.40530x*3 - 5.9470e-2x%4 


Surface Temperature 178 
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Surface Temperature (deg C) 


200 


100 


y = 180.95 + 20.243x - 0.66705x*2 - 0.19116x%3 - 3.9773e-2x%4 


Surface Temperature 181 
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Surface Temperature (deg C) 


y = 192.90 + 17.843x - 2.2433x%2 - 0.57151x*3 - 6.1657e-2x%*4 


200 


Surface Temperature 193 


100 
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Surface Temperature (deg C) 


y = 207.76 + 17.544x - 2.0301x*2 - 0.50866x%3 - 5.0198e-2x%4 


300 


Surface Temperature 208 


200 


100 
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Quartz: 


Thermal Conductivity: 


Specific Heat: 


Density: 


Macor: 


Thermal Conductivity: 


Specific Heat: 


Density: 


Aluminum: 


Thermal Conductivity: 


Specific Heat: 


Density: 


Material Properties: 


1.758 W/m-K 
812.2 J/kg-K 


2202 kg/m? 


1.297 W/m-K 
888.9 J/kg-K 


2520 kg/m? 


121.17 W/m-K 


962.32 J/kg-K 


2771 kg/m? 
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QuaLlezesolra 


EVAPORATION TIMES (sec) 


Volume 50 40 30 20 10 
Surface Temp (C) 
94 22514 205.8 auis}P4g ah WIA) 6 @ 88 
103 TS 33 163.9 13 Say7 OWA. f 70.6 
ita! 164.7 141.6 elviete 95 6779355 
124 3 116.9 97.6 of 525 
144 99.5 87 69.7 54.5 CHES 9/ 
139 102.889 89 79.5714 61-5 44.5 
158 82.4 USS 7 62.1 47.9 3 2 ter 
165 79.08 69.09 58.4 46.9 Shaver 
179 61 56 45.5 393 PAG 22 
190 Si) 49.5 Lada Fi BWV a PAB) 5 F) 
220 Vac Fi 36.7 3372 5 ak Zi > 18.8 
240 33.4: 30.82 8022 Pipks al 14 
154 89.45 76.64 65.2 S05 S38! 3877/3 
SNS 833 oes Osi: 48.54 3463 
142 99.8 92.6 eth 58.6 BS 
LYS, 102.9 87.8 80.5 66.7 44.5 
doai2 OW) 89.4 Poles 69.2 46.5 
128 114 Ones, S23 66.7 SOne 
119 243 0: alps s} COZ, 71.4 Sais 
98 PLOT As 188.4 156.5 WAS a 97.8 
BETA ZERO 
Surface Temp (C) 50 40 30 20 10 AVG 
94 1S 607) 1.559 P53 5 5/5 ab 1.5482 
103 ak 5 Syste) L529 eS Die abn Sys} Bhs SIS YS! 1.5346 
abals| 1.564 ahs feyeyal 1.548 She sSyoy:} ib a doyalak 1.5464 
124 aL slang es: 1.526 Heeslajsjal Gg SEY eos 1.542 
139 hs Syz! 1g BI} abo isis) ab BIS} 1.56 1542 
144 1.504 1.486 e547 die 539 505 tpl 
158 a Syy/ 524 1 oS 1.493 AL era? aka ial s3) 
165 UGS 1.496 1.492 1.495 thes S33) V5026 
154 he EYE bales AN, 5 $35 156 ib Gaal. 
53 Ibo Shs: 1.698 O57, 1.607 1.43 
142 1.594 a aCye 1.636 ib, Sys) 55 
LsiZ 1.64 1.66 1.68 1.566 ep o 2 
128 Ua ES 1.625 1.6 1.599 1.599 
1i9 1.66 bis 7 1.65 1.68 1.66 
98 1.68 1.62 1.596 1.66 1.64 
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Additional Macor Tests (Infrared Thermography Runs) 
Evaporation Time (sec) 


50 30 10 
Surf Temp (C) 

fest Po 146 73 
124 124 99 46 
143 94 wha 56 
SZ 82 62 30 
160 70 50 28 
178 46 28 16 
181 42 25 16 
Nie) 28 20 14 
208 Zs digi it 


BOILING REGIME EVAPORATION TIME (SEC) 


QUARTZ ; MACOR 


Initial Solid Surface Temperature ( C) 
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OO MNONSZWHND 


—~ 
_ 


<< <8 adi at ot ot 
NOW WP 


Surface Temp 
Macor 
160.000' 


155.000 
140.000 
133.000 
119.000 
115.000 
10S$.000 
102.000 
96.000 
Aluminun 
100.000 
95.000 
91.000 
87.000 
82.009 
78.000 
75.000 


On OWeB WR — 


RAYLEIGH # 


82700.000 
82290.000 
80460.000 
79170.000 
75600.000 
74300.000 
72180.000 
69200.000 
66730.000 


37150.000 
35900.000 
34810.000 
33660.000 
32090.000 
30730.000 
29700.000 


ime 


0.000 
1.000 
10.000 
20.000 
30.000 
40.000 
50.000 
60.000 
70.000 
80.000 
90.000 
100.000 
110.000 
120.000 
130.000 


a=.59 


12.270 
12.180 
11.890 
11.720 
11.400 
11.300 
11.130 
10.920 
10.730 


9.410 
9.270 
9.160 
9.030 
8.870 
8.730 
8.620 


oR 


124.000 


115.000 
118.000 
120.000 


a=.84 


17.470 
17.340 
16.920 
16.690 
16.230 
16.080 
15.840 
15.540 
15.280 


13.400 
13.210 
13.040 
12.860 
12.630 
12.420 
12.260 


1R 


124.000 
44.000 
73.800 
75.000 
78.600 
78.600 
77.400 
79.800 
81.000 
86.700 
89.000 
99.900 
116.000 
118.000 
120.000 
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a=.54 


135230 
11.150 
10.8560 
10.730 
10.430 
10.340 
10.180 

9.990 

9.820 


8.620 
8.490 
8.380 
8.270 
8.120 
7.990 
7.880 


transient data 


HC DATA 
a=.72 az.70 
14.970 14.560 
14.860 14.450 
14.510 14.110 
14.310 13.910 
13.910 13.530 
13.790 13.400 
13.580 13.200 
13.320 12.950 
13.100 12.730 
11.490 W1a170 
11.320 11.000 
WiAsNZ70 10.860 
11.020 10.720 
10.820 10.520 
10.650 10.350 
10.510 10.220 
2R 3R 
124.000 124.000 
120.000 124.000 
111.000 119.000 
114.000 118.000 
110.000 117.000 
111.000 118.000 
113.000 117.000 
112.000 118.000 
108.000 117.000 
108.000 116.000 
108.000 116.000 
113.000 117.000 
116.000 118.000 
118.000 119.000 
120.000 120.000 


4R 


124.000 
124.000 
122.000 
122.000 
120.000 
121.000 
119.000 
120.000 
120.000 
120.000 
119.000 
119.000 
119.000 
120.000 
123.000 


Experimental 


12.930 
12.000 
11.000 
19.250 
11.480 
11.540 
11.000 
11.000 
11.450 


11.800 
10.580 
10.670 
10.140 
10.420 

9.760 

8.710 


5R 


124.000 
124.000 
124.000 
124.000 
121.000 
121.000 
121.000 
121.000 
121.000 
121.000 
121.000 
121.000 
121.000 
123.000 
124.000 


Radius of Influence data 


Droplet Volume 


50 30 10 Evaporation time Surface Temp (C) 
SS 78 160 
7.79 67 160 

(ahpyip Ws 53 160 

8.74 Sy 160 

4.48 27 160 

5.8 30 160 

9.78 3h7/ 143 

Was 36 143 

Ws eis 75 143 

6.79 74 143 

Gaol 96 143 
6.09 93 143 
Gye a! 120 124 
Gree 128 124 
ere. 97 124 

ies OS 101 124 

6.97 73 101 

7) a PAS 139 101 

slegsis| 150 101 

6.06 195 101 
6.96 182 101 
LORS 48 24 

9.8 49 124 

6.9 82 aba 
6.97 85 LS 2 
Ow 63 L52 

USF 65 Siz 

Bices 3. auisy72 

Sieai2 26 178 
Se 29 178 
4.76 18 178 

Fhiesisya| ay) 208 
Lwee 19 193 
2.98 23 193 
JO2 14 193 

6.6 16 181 

(yo 25] 16 181 

Phe Gh 25 181 
5.66 24 18] 
5 30 181 
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